Duality considerations about the Maxwell-Podolsky theory through 
the symplectic embedding formahsm and spectrum analysis 
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We find the dual equivalent (gauge invariant) version of the Maxwell theory in D = 4 with 
a Proca-like mass term by using the symplectic embedding method. The dual theory obtained 
(Maxwell-Podolsky) includes a higher-order derivative term and preserve the gauge symmetry. We 
also furnish an investigation of the pole structure of the vector propagator by the residue matrix 
which considers the eventual existence of the negative-norm of the theory. 
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I. INTRODUCTION 

The interest in the dual mapping between two different 
theories that show the same physical properties has been 
increased in the last years motivated by its success in 
both field theory as well as string theory At the 
same time, the interest in he study of theories involving 
higher-order derivatives is by now well appreciated and 
remains intense. Within the context of Maxwell theory, 
generalizations involving higher-order derivatives can be 
found in dig. 

We know that in _D — 3 + 1 massive gauge fields are 
introduced by the Higgs Mechanism through the sponta- 
neous breaking of symmetry. However, in the Abelian 
case, it is possible to introduce a massive gauge field 
without having to accomplish a spontaneous breaking of 
symmetry. So, the gauge symmetry does not forbid the 
appearance of the massive gauge fields. We can fall back 
to a topological term of mass like Chern-Simons Q, that 
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in D ~ 2 + 1 preserve the gauge and Lorentz symmetries. 
However, in Z) = 3 -I- 1 the Lorentz symmetry is broken 
for the Chern-Simons theory. 

There are different ways of extending electrodynamics, 
trying to smooth infrared or ultraviolet singularities that 
appear at large or short distances. One is the Born-Infeld 
2] type of generalization, that involves nonlinear exten- 
sion. The other includes the generalizations introduced 
by Proca that involves the addition of a mass term 
for the vector field, which was introduced to smooth in- 
frared singularities. But, the gauge symmetry is clearly 
lost, due to the term A^A'^. Another way is through the 
work of Podolsky Q , that involves higher-order deriva- 
tives and was introduced to smooth ultraviolet singular- 
ities. The Podolsky theory propagates a massive photon 
and the Lorentz and gauge symmetries are not sponta- 
neously broken. Besides, it has a fundamental role in 
the discussion about the compatibility of the magnetic 
monopole and massive photons. 

Working in D = 3, the authors in [sj presented 
the dual mapping of the Maxwell-Chern-Simons-Proca 
(MCS-Proca) with a MCS-Podolsky term through the 
iterative gauge embedding procedure It was intro- 
duced a convenient parameter a in the Maxwell term. 
When a = we have a MCS-Proca, that is the stan- 
dard self-dual model , and when a = 1 it was obtained 
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an extended CS-Podolsky model. The objective of the 
introduction of this parameter is to construct a master 
action, from which one can obtain both the model and 
its dual partner. But one can question if this dual part- 
ner is unique. We believe that the symplectic embedding 
technique can furnish a whole family of dual models be- 
cause of the arbitrariness of the zero-mode which is also 
the generator of the infinitesimal gauge transformations, 
as we will see below. 

Consequently, the purpose of this paper is to work out 
the dual aspects of electrodynamics that appears in the 
Proca and Podolsky models by the symplectic embedding 
method @, that has been very efficient to reveal a de- 
sired hidden symmetry and a dual partner of the model 
0) [13, [IHj We analyze the residue matrix at each 
pole of the propagators. The poles furnishes a relation 
between energy and momentum which can be associated 
to a massive particle. The residue matrix gives informa- 
tion about the degree of freedom of the partner polariza- 
tion, where physical states are partnered with positive 
eigenvalues of the residue matrix at each pole [1^. The 
existence of negative eigenvalues shows that in these sit- 
uations nonphysical states that correspond to negative 
norm particle states (ghosts) are introduced. 

The symplectic embedding method (SEM) ^] is not 
affected by ambiguity problems. It has the great advan- 
tage of being a simple and direct way of choosing the 
infinitesimal gauge generators of the built gauge theory. 
As a consequence, we have the freedom to choose the con- 
tent of the embedded symmetry according to our neces- 
sities. This feature makes possible a greater control over 
the final Lagrangian. This method can avoid the intro- 
duction of infinite terms in the Hamiltonian of embedded 
non-commutative and non-Abelian theories. This can be 
accomplished because the infinitesimal gauge generators 
does not result from previous unclear choices. Another 
advantage related is the possibility of doing a kind of 
general embedding. In other words we can say that, in- 
stead of choosing the gauge generators at the beginning, 
one can leave some unfixed parameters with the aim of 
fixing them later, when the final Lagrangian has being 
achieved. Although a faster way to obtain the final the- 
ory is through the fixing of such parameters, it is more 
interesting in order to obtain knowledge of the in ques- 
tion. Also, it is helpful if we want to know its hidden 
symmetries, but some aspects of the Lagrangian are re- 
quired. The embedding approach is not dependent on 
any undetermined constraint structure and it also works 
for unconstrained systems. This is different from all the 
existent embedding techniques that we can use to con- 
vert [13, [lB|, to project [l6| or to reorder [13] the exis- 
tence of the second-class constraints into a first-class sys- 
tem. This technique on the other hand only deals with 
the symplectic structure of the theory, so this embedding 
structure does not rely on any pre-existent constrained 
structure. 

In order to make this paper self-sustained, in the next 
section we present a brief review of the dual embedding 



formalism. In the section IIIIl the SEM will be used in 
the Proca model to construct a massive gauge invariant 
theory. In the section IIVI we present an investigation 
of the structure of the gauge propagator. Finally, the 
conclusions are accomplished in the last section, as usual. 

Our metric tensor has signature (H ) and we use 

natural units. 



II. THE DUAL EMBEDDING FORMALISM 

This technique follows the Faddeev-Shatashivilli's sug- 
gestion jl8j and is set up on a contemporary framework 
to handle constrained models, i.e., the symplectic for- 
malism [l9|. In the following lines, we try to keep this 
paper self-sustained reviewing the main steps of the dual 
embedding formalism. We will follow closely the ideas 
contained in [Tl| . 

Let us consider a general noninvariant mechanical 
model whose dynamics is governed by a Lagrangian 
C{ai, cLi, t), (with 1 = 1,2,..., TV), where and di are the 
space and velocities variables, respectively. Notice that 
this model does not result in the loss of generality nor 
physical content. Following the symplectic method the 
zeroth-iterative first-order Lagrangian one-form is writ- 
ten as 



and the symplectic variables are 



(0)6 ^ ^ 

Pi 



with 6 
with 6 



1,2,. ..,N 
N + l,N + 2,. 



,2N, 



(1) 



(2) 



where A'"g ^ are the canonical momenta and V^^'' is the 
symplectic potential. From the Euler-Lagrange equations 
of motion, the symplectic tensor is obtained as 



^(0) 
J 0/3 



(0) 



(3) 



If the two-form / = ^fepd^^ A d^^ is singular, the sym- 
plectic matrix ([3]) has a zero- mode (i/(°)) that generates 
a new constraint when contracted with the gradient of 
the symplectic potential. 

This constraint is introduced into the zeroth-iterative La- 
grangian one-form equation ((T|) through a Lagrange mul- 
tiplier 77, generating the next one 

= AWde^i^^-F(i)(C)di, (5) 
with 7 = 1,2,..., (2iV + 1) and 

= (^W^,7), (6) 

41) = {4'\n("^ 
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As a consequence, the first-iterative symplectic tensor is 
computed as 



Ai) _ 

■1713 - 



dA. 



(1) 



(7) 



If tliis tensor is nonsingular, the iterative process stops 
and the Dirac's brackets among the phase space variables 
are obtained from the inverse matrix (/^Jj*)"^ and, con- 
sequently, the Hamilton equation of motion can be com- 
puted and solved, as discussed in [20]. It is well known 
that a physical system can be described at least clas- 
sically in terms of a symplectic manifold M . From a 
physical point of view, M is the phase space of the sys- 
tem while a nondegenerate closed 2-form / can be iden- 
tified as being the Poisson bracket. The dynamics of the 
system is determined just specifying a real-valued func- 
tion (Hamiltonian) H on the phase space, i.e., one of 
these real- valued function solves the Hamilton equation, 
namely, 



(8) 



and the classical dynamical trajectories of the system in 
the phase space are obtained. It is important to mention 
that if / is nondegenerate, equation has an unique 
solution. The nondcgeneracy of / means that the lin- 
ear map b : TM T*M defined by \){X) := \?{X)f is 
an isomorphism, due to this, the equation ([5]) is solved 
uniquely for any Hamiltonian (X = [>~^{dH)). On the 
contrary, the tensor has a zero-mode and a new constraint 
arises, indicating that the iterative process goes on un- 
til the symplectic matrix becomes nonsingular or singu- 
lar. If this matrix is nonsingular, the Dirac's brackets 
will be determined. In Ref. [20,], the authors consider 
in detail the case when / is degenerate. The main idea 
of this embedding formalism is to introduce extra fields 
into the model in order to obstruct the solutions of the 
Hamiltonian equations of motion. We introduce two ar- 
bitrary functions which are dependent on the original 
phase space and of WZ's variables, namely, 5'(ai,pi) and 
G{ai,pi,r]), into the first-order Lagrangian one- form as 
follows 



with 



F(o) +G(a„p„r7), (10) 



where the arbitrary function G{ai,pi,ri) is expressed as 
an expansion in terms of the WZ field, given by 

oo 

G(a„p„7/)-^g(")(a„p„r/), g(")(a„p„ r;) ^ 77" , 

n=l 

(11) 

and satisfies the following boundary condition 

G(a„p„7y = 0) = 0. (12) 



The symplectic variables were extended to also contain 
the WZ variable = (C^°)^ ry) (with ^=1,2,..., 2N+ 
1) and the first-iterative symplectic potential becomes 

00 

n=l 

In this context, the new canonical momenta are 
7(0) _ / with 0~ =1,2,. .. ,2N 



^f, with 61== 2N -I- 1 



and the new symplectic tensor, given by 



(14) 



■'9/3 



dA 



(0) 



a|(o)e a|(o)/3 ' 



that is 



/AO) AO)^ 
f(o) _ / Jef3 Jen 
■lep ~ AO) 

\Jr,l3 







(15) 



(16) 



To sum up we have two steps: the first one is addressed 
to compute 5'(ai,pi) while the second one is dedicated to 
the calculation of G{ai,Pi,ri). In order to begin with 
the first step, we impose that this new symplectic ten- 
sor (/^°^) has a zero-mode i>, consequently, we get the 
following condition 



r,{0)6A0) ^ 
•I or-, 



0. 



(17) 



At this point, / becomes degenerate and, in consequence, 
we introduce an obstruction to solve, in an unique way, 
the Hamilton equation of motion given in equation ([8]). 
Assuming that the zero-mode 



K 1) , 



(18) 



and using the relation given in (fT7| together with 
we get a set of equations, namely. 



where 



J) AO) AO) _ Q 



,(0) _ dAf 



/,,/3 - 



(19) 



(20) 



The matrix elements fi^ are chosen in order to disclose 
a desired gauge symmetry. Note that in this formalism 
the zero-mode z/C^^ is the gauge symmetry generator. At 
this point, it is worth to mention that this characteristic 
is important because it opens up the possibility to dis- 
close the desired hidden gauge symmetry from the non- 
invariant model. It awards to the symplectic embedding 
formalism some power to deal with noninvariant systems. 
From relation (|17p some differential equations involving 



4 



^{ai,pi) are obtained, equation (fT9l) . and after a straight- 
forward computation, "^{ai^Pi) can be determined. 

In order to compute G{ai,pi, rf) in the second step, we 
impose that no more constraints arise from the contrac- 



tion of the zero- mode {v^'^^^) with the gradient of the po- 
tential V^^\a^,p^ , rf). This condition generates a general 
differential equation, which reads as 



~{0)e dV^^\auP„r^) 
^ 9|(o)e 



977 



drj 



0, 



... = 



(21) 



(22) 



that allows us to compute all correction terms 
Q'-'^\ai,pi,r]) in order of r]. Note that this polynomial 
expansion in terms of rj is equal to zero, consequently, 
whole coefficients for each order in 77 must be null iden- 
tically. In view of this, each correction term in order of 
77 is determined. For a linear correction term, we have 



o dV(°Ha,,p,) dg(^Ha,,p,,rj) _ 

For a quadratic correction term, we get 

,dg('Ha^,p,,il) , dg'^^\a,,p,,7j) 



(23) 



0. (24) 



From these equations, a recursive equation for 77 > 2 is 
proposed as 



e dg("-^Ha,,p,,r^) , 9g(")(a„p„ 77) _^ 



drj 



(25) 



that allows us to compute the remaining correction terms 
in order of 77. This iterative process is successively re- 
peated until (pij) becomes identically null, consequently, 
the extra term G{ai,pi,ri) is obtained explicitly. Then, 
the gauge invariant Hamiltonian, identified as being the 
symplectic potential, is obtained as 

H(a„p„77) = F(°)(a„p,) + G(a„p„ 77), (26) 

and the zero-mode is identified as being the gener- 
ator of an infinitesimal gauge transformation, given by 



5e' = e^^('"', 
where e is an infinitesimal parameter. 



(27) 



III. MAXWELL-PODOLSKY FROM 
MAXWELL-PROCA 

In this section we follow the steps described in the 
last one in order to find a dual equivalent action to the 



original theory. As we will see this dualization allows us 
to construct a mapping between a non gauge invariant 
model (Proca) and a gauge invariant one (Podolsky). 

Let us consider the massive Lagrangian density in four 
dimensions 



/:MProca = -Ji^^.F''" + |A^A^, (28) 

where the canonical momenta are 

TT, = -Fq,. (29) 

Hence, we have the following Lagrangian, 

£(0) = n'A, ~ (30) 

where 

= ^F,,F^^ - ^A^A^^ - i^V, + tt.SMo. (31) 

Considering a convenient zero-mode as, 

i>(°) = ( a" -1 ) (32) 

and extending the phase space with the introduction 
of the WZ fields and the two arbitrary function ifj = 
{A\ ij\ A°,r]) and G = G {A\ ij\ A° , rj) we change the 
Lagrangian ([50]) . namely. 



£(0) = n'A, +^Pr]- 
where the symplectic potential is 



(33) 



= ^F,,F'^ - ^A^A^^ ~ IttV, + tt.SMo + G. (34) 
The function G is expressed as 

oo 

G{A\^\A'',7]) ^ith g"cx77". (35) 

n=0 
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Now, contracting the zero-mode (p2)) with the symplec- 
tic matrix of the Lagrangian and after a straight- 

forward computation, we get 

^ = -dy. (36) 

Then, the Lagrangian becomes 

=^*i, -a.vr^Ty-y^"). (37) 

Contracting the zero-mode ((32|) with the symplectic 
potential we have that, 

j dy{[~dlF,,{y) ~ ^A,{y)] a^<5(x - y) 

+ [-a>.(y) - My)] d"j{^ - y) - ^^^) = 

(38) 

with the solution 
0(0) = {-d'F.j - ^Aj) d^rj + (-9V, - e^o) a°ry. (39) 
With this result we have a new symplectic potential, 

oo 

yd) = + ^(0) + J2 (40) 

and again, a contraction of this result with the zero-mode 
give us the next correction to the function G: 

J dy {-^dlr,iy)d^S{^ - y) 

- ^d^jViy)dSS{^ - y) - ^^^) = 0, (41) 

namely, 

gW ^ -^d.Tjd^Tj - ^dorjd^rj. (42) 

Note that the second-order correction term has depen- 
dence only on the WZ field, thus all the correction terms 
^;(") for n > 2 are zero. Then, the gauge invariant first- 
order Lagrangian density, after some algebra, is given by 

C = CMProca + J^d^V + fs^^S^T? (43) 

where CMProca is given in ((28)) and 

= d-'F,^ + M (44) 

is the Euler current. 

For convenience, let us define 9^77 as an external field 
Bfj , and we can rewrite (|^ as 

C = CMProca + Jt^B'' + ^Bf,B^' , (45) 



and solving the equation of motion for B^ we have, 

C = CMProca — 2^"^A'"^^- (^6) 

Finally, eliminating the WZ fields we can work alge- 
braically to obtain the final equivalent theory as 

^ = -\f,^F'^'' ~ Y^idxFnid'F,^), (47) 

which is the Maxwell-Podolsky model [if], where the di- 
mension of ^ is [mass]'^. 

We found that, starting from the Maxwell-Proca 
model, where the gauge symmetry is not present, and 
using the symplectic embedding method we found the 
Maxwell-Podolsky model that propagates a massive pho- 
ton and preserves the gauge symmetry. 

As well known from the symplectic formalism litera- 
ture, the zero-mode is the generator of the infinitesimal 
gauge transformation given by SO — ei>(°). So, using 
the zero-mode given in (I32p . we have the following gauge 
symmetries, 

SA, = -a,e 
Stt, = 
Mo = ~do€ 
Srj = e, (48) 

where e is an infinitesimal time-dependent parameter. 

The gauge transformations obtained above completes 
the process of symplectic embedding which is altogether 
different from the gauge embedding technique, as stressed 
above. We believe that the arbitrariness of the zero-mode 
choice in order to obtain a whole family of dual partners 
is a positive factor when we compare both procedures. 

It is important to notice the dual aspects of the electro- 
dynamics that appear in the Proca and Podolsky mod- 
els. This duality is also important in the extension of the 
bosonization programme from _D 2 to higher dimen- 
sions [1]. 

IV. THE SPECTRUM ANALYSIS OF THE 
MAXWELL-PODOLSKY THEORY 

In this section, we derive the propagators of the 
Maxwell-Podolsky theory in order to carry out a study 
about its spectrum. 

After a straightforward computation, we can write the 
Lagrangian (liT)) as 

(49) 

Let us rewrite (|49|) conveniently as, 

C = iA^e^M^ (50) 
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where 

= - df,d^ ~ ^g^,^a'^ + ^adf,d^ (51) 

is a differential operator to the theory with g^^ = 6^^ + 
(jj^i,. The is the transversal projector and cj^i, = 
is the longitudinal projector. 

The inverse operator can then be written as 



where the vector propagator is defined as 



(52) 



(53) 



In momentum space the (1551) is expressed by the relation 



{AM 



{AM 



—ft 



11.2 i^^^ 



(54) 



(55) 



The relation (|55p is the vector propagator for the Maxwell 
theory extended due to the Podolsky term in the (|47)) . 
Analyzing the dispersion relations corresponding to the 
propagator of the Maxwell-Podolsky, ([T71) . we find one 
particle of spin 1 (one massless pole = 0) and other 



particle of spin 1 with a massive pole k 



If e < 



then fc^ = jiy is a massive "photon" and it is not a 
tachyon. 

Now, let us to analyze the eigenvalue of the residue 
matrix for the pole k^ — —j with ^ < , choosing /c^ 

purely timelike /c^ = (|^| ' , 0), we get 

k,, k,. 



k^ku 



or 



R = 



/O ON 

i 

i 

VO 



(56) 



(57) 



Hence, for timelike fc'', the poles of {A^j^A^) respect 
causality (they are not tachyonic) and correspond to 
physically acceptable one-particle states with three de- 
gree of freedom, since the residue matrix exhibits a sin- 
gle positive eigenvalue. Consequently, the model may be 
adopted as a consistent theory. 

V. FINAL DISCUSSIONS 

In this work we used the symplectic embedding pro- 
cedure to promote the dualization of the gauge invariant 
Maxwell theory in 13 = 4, modified by the introduction of 
an explicit massive (Proca) term. This technique was in- 
troduced originally to deal with constrained models and 
involves the construction of a symplectic tensor and a 
zero-mode term. 

The new theory obtained (the dual partner) involves a 
massive term for the vector field, which is the Podolsky 
term that involves higher-order derivatives and propa- 
gates a massive "photon" without a spontaneous break- 
ing of the gauge symmetry. 

In comparison with the gauge embedding formalism, 
although being a swift iterative technique in order to 
obtain a dual partner of the original theory, a question 
about the uniqueness lingers. To use the symplectic em- 
bedding technique means to obtain a whole family of dual 
partners thanks to the arbitrariness of the zero-mode 
term. Besides, the zero-mode term is the generator of 
the infinitesimal gauge transformations of the dual part- 
ner theory. 

The duality between the Proca and the Podolsky theo- 
ries is also important in the extension of the bosonization 
programme from D = 2 to higher dimensions. Although 
out of the scope of this work, we can say that bosoniza- 
tion is very important for the interpretation of the new 
parameters in the Podolsky extension. 

To infer about the physical nature of the simple poles, 
we have to calculate the eigenvalues of the residue matrix 
for each of these poles. To carry out this investigation 
we analyzed the pole structure of the vector propaga- 
tor by the residue matrix, where we found three positive 
eigenvalues that correspond to the physically acceptable 
one-particle states with three degree of freedom. When 
^ — + 0, we easily obtain the original QED (massless pho- 
ton). 
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